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Abstract 

We study the interactions between the /o(980) and ao(980) scalar res- 
onances and the lightest pseudoscalar mesons. We first obtain the 
elementary interaction amplitudes, or interacting kernels, without in- 
cluding any ad hoc free parameter. This is achieved by using previous 
results on the nature of the lightest scalar resonances as dynamically 
generated from the rescattering of S-wave two-meson pairs. After- 
wards, the interaction kernels are unitarized and the final S-wave am- 
plitudes result. We find that these interactions are very rich and 
generate a large amount of pseudoscalar resonances that could be as- 
sociated with the #(1460), vr(1300), vr(1800), 7/(1475) and Jf (1835). 
We also consider the exotic channels with isospin 3/2 and 1, having 
the latter positive G-parity. The former could be also resonant in 
agreement with a previous prediction. 
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1 Introduction 



Due to the spontaneous chiral symmetry breaking of strong interactions [IHI] 
strong constraints among the interactions between the lightest pseudoscalars 
arise, which are most efficiently derived in the framework of Chiral Pertur- 
bation Theory (CHPT) [SHE]. For the isospin (J) 0, 1 and 1/2 the scattering 
of the pseudoscalars in S-wave is strong enough to generate dynamically the 
lightest scalar resonances, namely, the /o(980), a (980), k and a, as shown 
in refs. [9T4T3] . Still one can make use of the tightly constrained interactions 
among the lightest pseudoscalars in order to work out approximately the scat- 
tering between the latter mesons and scalar resonances, as we show below. 
We concentrate here on the much narrower resonances /o(980) and ao(980) 
and consider their interactions with the pseudoscalars ir, K, 77 and rf '. If these 
interactions are strong enough new pseudoscalar resonances with J PC = 0~ + 
would come up. This is the case and the resulting pseudoscalar resonances 
have a mass larger than 1 GeV (this energy limit is close to the masses of 
the / (980) or a (980)), typically following the relevant scalar-pseudoscalar 
thresholds. 

The problem of the excited pseudoscalars above 1 GeV is interesting by 
itself. These resonances are not typically well-known [13]. In / = 1/2 one 
has the K(14Q0) and 7^(1630) resonances. The 1 = 1 resonances 7r(1300), 
7r(1800) are somewhat better known [13]. They are broad resonances with a 
large uncertainty in the width of the former, which is reported to range be- 
tween 200-600 MeV in the PDG [13] • Some controversy exists for interpreting 
the decay channels of the 7r(1800) within a quarkonium picture [T5lfT6] . It was 
suggested in [T5] that together with the second radial excitation of the pion 
there would be a hybrid resonance somewhat higher in mass [751 [TT] . Spe- 
cial mention deserves the 1 = channel where the 7/(1295), 7/(1405), 7/(1475) 
have been object of an intense theoretical and experimental study. For an 
exhaustive review on the experiments performed on these resonances and the 
nearby 1 ++ axial- vector resonance /i(1420) see ref. [18J. Experimentally it 
has been established that while the 77(1405) decays mainly to ao7r the 7/(1475) 
decays to K*K + c.c [THUS]. m this way, the study of the rj-n-n system is 
certainly the most adequate one for isolating the 7/(1405) resonance because 
both the /i(1420) and 7/(1475) have a suppressed partial decay width to this 
channel [13]. Refs. [T4T[T8] favor the interpretation of considering the 7/(1295) 
and 7/(1475) as ideally mixed states (because the 7/(1295) and the 7r(1300) are 
close in mass) of the same nonet of pseudoscalar resonances with the other 
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members being the 7r(1300) and X(1460). All these resonances would be the 
first radial excitation of the lightest pseudoscalars [TJ]. The 77 (1405) would 
then be an extra state in this classification whose clear signal in gluon-rich 
process, like pp [191120] or J/ty radiative decays [2T1[22], and its absence in 
77 collision [23J, would favor its interpretation as a glueball in QCD [241125]. 
However, this interpretation opens in turn a serious problem because present 
results from lattice QCD predict the lowest mass for the pseudoscalar glueball 
at around 2.4 GeV [26T28J. Given the success of the lattice QCD prediction 
for the lightest scalar glueball, with a mass at around 1.7 GeV [29l[30], this 
discrepancy for the pseudoscalar channel would be quite exciting. QCD sum 
rules [31] give a mass for the lightest pseudoscalar gluonium of 2.05±0.19 GeV 
and an upper bound of 2.34 ±0.42 GeV. However, the 77(1405) would fit as a 
0~ + glueball if the latter is a closed gluonic fluxtube [32]. On the other hand, 
it has also been pointed out that the mass and properties of the 7/(1405) are 
consistent with predictions for a gluino-gluino bound state [25j[331|3I]. The 
previous whole picture for classifying the lightest pseudoscalar resonances has 
been challenged in ref. [16] . The authors question the existence of the 7/(1295) 
and argue that, due to a node in the 3 P wave function of the //(1475) [35J, 
only one isoscalar pseudoscalar resonance in the 1.4-1.5 GeV region exists. 
This node shifts the resonant peak position depending on the channel, aovr 
or K*K + a a It has been also recently observed by the BES Collaboration 
the resonance X(1835) with quantum numbers favored as a pseudoscalar 0~ + 
resonance both in J/\I> — » ^ypp |36j and in J/$ -)■ 77r + 7r~7/ [37]. For the for- 
mer decay ref. [38J offers an alternative explanation in terms of the pp final 
state interactions. 

We consider here the S-wave interactions between the scalar resonances 
/o(980) and a (980) with the pseudoscalar mesons tc, K, 77 and r( . The 
approach followed is an extended version of the one that refs. [391140] applied 
to study the S-wave interactions of the 0(1020) with the /o(980) and ao(980) 
resonances, respectively. We show that the interactions derived generate 
resonances dynamically that can be associated with many of the previous 
pseudoscalar resonances listed above, namely, with the K(1460), 7r(1300), 
7r(1800), 7/(1475) and X(1835). In this way, new contributions to the physical 
resonant signals result from this novel mechanism not explored so far. In 
addition, we also study other exotic channels and find that the 7 = 3/2 a^K 
channel could also be resonant. 

After this introduction we present the formalism and derive the S-wave 
scattering amplitudes for scalar-pseudoscalar interactions in section [2] Sec- 
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Figure 1: Triangle loop for calculating the interacting kernel for 
Si(pi)P\(ki) — > 52(^2)^2(^2), where between brackets the four-momentum 
for each particle is given. Sip represents the initial, final scalar resonances 
and similarly for regarding the pseudoscalar mesons. 

tion [3] is dedicated to present and discuss the results. Conclusions are given 
in section HI 

2 Formalism. Setting the model 

Our approach is based on the triangle diagram shown in Fig. [T] where an 
incident scalar resonance S\ decays into a virtual KK pair. The filled dot 
in the vertex on the bottom of the diagram corresponds to the interaction of 
the incident (anti)kaon in the loop with the pseudoscalar Pi giving rise to the 
pseudoscalar P 2 and the same (anti)kaon. The out-going scalar resonance is 
denoted by S2. The basic point is that this diagram is enhanced because 
the masses of both the /o(980) and ao(980) resonances are very close to 
the KK threshold. In this way, for scattering near the threshold of the 
reaction, one of the kaon lines in the bottom of the diagram is almost on- 
shell. Indeed, at threshold and in the limit for the mass of the scalar equal 
to twice the kaon mass this diagram becomes infinite. This fact is discussed 
in detail in ref. [39] where it was already applied for studying successfully 
the 0(lO2O)/o(98O) scattering and the associated 1 F(2175) resonance. 
The BABAR [H] and BELLE 02] data on e+e" </>(1020)/ (980) were 
reproduced accurately, where a strong peak for the latter resonance arises. 
An important conclusion of [39J is that the Y (2175) can be qualified as being a 
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resonance dynamically generated due to the interactions between the 0(1020) 
and the / (980) resonances, see also ref. [43J. This work was extended to 
/ = 1 in jlQ] for studying the 0(lO2O)ao(98O) S-wave. There it was remarked 
the interest of measuring the cross sections e + e~ — > 0(lO2O)7r°?7 because it 
is quite likely that an isovector companion of the y(2175) appears. In our 
present study, as well as in refs. [321110], one takes advantage of the fact 
that both the / (980) and a (980) resonances are dynamically generated by 
the meson-meson self-interactions J9j[13j|ll]. This conclusion is also shared 
with other approaches like refs. [13116]. In this way, we can calculate the 
couplings of the scalar resonances considered to two pseudoscalars, including 
their relative phase. The coupling of the / (980) and a (980) resonances to 
a KK pair in / = and 1, respectively, is denoted by gj and g ao . These 
states \KK)i =Q and \KK)i=\ are given by 

\KK) I=0 = -^=\K + K-+K°K°) , 

\KK) I=l = — ^\K + K~ - K°K°) . (2.1) 
V2 

In this way, the / (980) couples to K + K-(K°K°) as --^(--^)g fo while the 

a (980) couples as -75(73)^0- 

Let us indicate by P the total four-momentum P = p\ + k\ = P2 + &2 in 
Fig. [TJ This diagram is given by g\giL>K, with g\ and g% the coupling of the 
initial and final scalar resonance to a KK pair, respectively, and Lk is given 
by 



dH T((P - 



(2tt) 4 {£ 2 -m 2 K + ie)(( Pl - £) 2 - m\ + ie)((p 2 - £) 2 - m 2 K + is) ' 

(2.2) 

In this equation T((P — £) 2 ) represents the interaction amplitude between the 
kaons with the external pseudoscalars. Here, we employ the meson-meson 
scattering amplitudes obtained in ref. [TJ] but now enlarged so that states 
with the pseudoscalar 7/ are included in the calculation of T((P — £) 2 ), as 
detailed in appendix [A] Interestingly, these amplitudes contain the poles 
corresponding to the scalar resonances a, k, / (980), a (980) and other poles 
in the region around 1.4 GeV [173]. 

In order to proceed further we have to know the dependence of T((P—£) 2 ) 
on its argument that includes the integration variable £. This can be done 



5 



by writing the dispersion relation satisfied by T(q 2 ) which is of the form 

rp/ 2\ rp/ x . (f ~ S A Resi q 2 -s A [°° , , ImT(s') 

T(g = T(s A ) + > -j + / ds — ^— . 

^ q 1 - Si Si- s A it J Sth {s' - q 2 ){s' - s A ) 

(2.3) 

One subtraction at sa has been taken because T(q 2 ) is bound by a constant 
for q 2 — > oo, with T(sa) the subtraction constant. Typically poles are 
also present deep in the g 2 -complex plane located at Sj whose residues are 
ReSj. These poles appear on the first Riemann sheet and are an artifact of 
the parameterization employed [T3"ll4"T] . For q 2 along the physical region they 
just give rise to soft extra contributions that could be mimic by a polynomial 
of low degree in q 2 . Inserting eq. (12. 3 j) into eq. (12. 2p . with (P — €) 2 = q 2 , one 
can write for 

L K = (t{s a ) + J2 J^O C 3 + V C,(s t )ReSi 
V s i — sa/ 



ds'lmT(s') 



C 3 



sa 



+ C 4 (s') 



(2.4) 



Here we have introduced the three- and four-point Green functions C 3 and 
C A (M 2 ) defined by 



C 3 = i 



dH 1 



(2tt) 4 (£ 2 -m 2 K + i£)((pi - £) 2 -m 2 K + ie){{p 2 - £) 2 - m 2 K + is) 

dH 1 

(2tt) 4 {£ 2 -m 2 K + ie)({ Pl - I) 2 - m\ + ie)((p 2 - £) 2 - m 2 K + ii) 

' (2-5) 



((P - £) 2 - M| + ze) 



Notice that M| can be real positive (when M| = s' in the dispersion relation), 
but it could also be negative or even complex when M| = Sj from the poles. 
One has still to perform the angular projection for C 3 and C^Mf). Once 
this is done, eq. (12 .4p can still be used but with C% and C±(Ml) projected 
in S-wave, as we take for granted in the following. These functions and 
their S-wave projection are discussed in appendix [Bj For Si{pi)Pi{ki) — > 
5*2(^2)^2(^2) we have the usual Mandelstam variables s = (pi + A4) 2 = (p 2 + 
h)\t = ( Pl -p 2 ) 2 = (h-k 2 ) 2 andu= ( Pl -k 2 ) 2 = (ps-Zd) 2 = M 2 Si + M 2 2 + 
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+ Mp 2 — s — t, with the masses of the particles indicated by M with the 
subscript distinguishing between them. The dependence on the relative angle 
9 enters in t as t = (p? - k^) 2 - (p-p') 2 = (p? - K) 2 -p 2 -p /2 + 2|p| |p'| cos6> 
with p and p' the CM three-momentum of the initial and final particles, 
respectively. 

Eq. (12. 4p is our basic equation for evaluating the interaction kernels. One 
has only to specify the pseudoscalars actually involved in the amplitude T(q 2 ) 
according to the specific reaction under consideration. We now list all the 
channels involved for the different quantum numbers and indicate the actual 
pseudoscalar-pseudoscalar amplitudes required as the argument of Lk- 

• i = o, G = +1 

2g 2 

T L (a 7c -> a 7c) = -^L K [AT^%l K - T'=^] , 

T L (a ir -»■ f rj) = 2gf gaoL K [T*K]ll K ] , 

T L (foV M = HMT^SIk] , 
T L (a 7r ->■ /oV) = 2g fo g ao L K [T I n r^ 7TK ] , 

Woi/ /oi/) = 2^^^^] . (2.6) 

• 7=1/2 



,2 



T L (f K -> /o-fC) = yLx[3T^ M + T k ^ kR \ , 
T L (foK a fO = y ^ 9 ^ 9ao L k [T k = r \ kR - T k =^ kR ] , 
T L (a K -> a ^) = ^^[3T^%^ + . (2.7) 
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• / = i, g = -1 



2o 2 

T L (f n -> /oTr) = -^^[2^-^ + T^ 2 ^] , 

T L (fo^a oV ) = 2 ^L K [T^ vK }, 

T L (a oV a 7?) = 2^ o L K [T;=^] , 
T L (/ovr^a V) = ^^[T^], 

Tl(o,q1] ->■ a V) = 2^ L^[T T ^,^] , 

T L (aoi/ -> <W) = [Tj?^] • (2.8) 



• / = 1, G = +1 

T L (a 7T aovr) = - t££k\ ■ (2-9) 



• 1 = 3/2 

Tl(o-qK a K) = 2^ ( L K [T^^] . (2.10) 



In the previous equations the different scalar-pseudoscalar states are pure 
isospin ones corresponding to the isospin 7 indicated for each item. This also 
applies to the pseudoscalar-pseudoscalar states, with I as indicated in the 
superscript of T. The symbol G refers to G-parity. On the other hand the 
1 = 3/2 ttK amplitude, being much smaller than the I = 1/2 one, has 
negligible effects, although it has been kept in the previous expressions. 




Figure 2: Iteration of the interaction kernels (denoted by the triangles in the 
figure) by inserting scalar(double lines)-pseudoscalar(dashed lines) interme- 
diate states. 



For each set of quantum numbers specified by the isospin I and G-parity 
G (if the latter is not defined this label should be omitted) we join in a 
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symmetric matrix Tig the different T^(i — > j) calculated above. Then, in 
order to resum the unitarity loops, as indicated in Fig. [2J and obtain the 
final S-wave scalar-pseudoscalar T-matrix, Tig, we make use of the equation 

Tig =[! + Tig • giois)}' 1 ■ Tig ■ (2.11) 



For a general derivation of this equation, based on the N/D method 
see refs. [131II9] and ref. [9], where it is connected with the Bethe-Salpeter 
equation. In eq. ( 12. lip gicis) is a diagonal matrix whose elements are the 
scalar unitarity loop function with a scalar-pseudoscalar intermediate state. 
For the calculation of giG(s)i, corresponding to the i t h state with the quantum 
numbers IG and made up by the scalar resonance Si and the pseudoscalar 
Pi, we make use of a once subtracted dispersion relation [13]. The result is 

If Mf. Mp_ - Mi + 8 , Ml 



giG{s)i = j^-r^ { a t + log —J- i — - — - log ■ 



+ 



(4tt) 2 [ 1 ' b fi 2 2s to Mp. 

log(s - A + 2Vs|p|) + log(s + A + 2Vs|p| 



IPI 



log(-s + A + 2Vs|p|) - log(-s - A + 2Vs|p|) 



(2.12) 



with |p| the three-momentum of the channel SiP{ for a given s and A = 
M Pi — Mg.. The subtraction at is restricted to have natural values so that 
the unitarity scale [39] ^irf^/ y/\a±\ becomes not too small (e.g. below the 
p-mass) so that |ai| < 3. In addition, we require the sign of at to be negative 
so that resonances could be generated when the interaction kernel is positive 
(attractive). 

As already indicated in ref. [10] to ensure a continuous limit to zero 
a (980) width, one has to evaluate Tig at the a (980) pole position with 
positive imaginary part so that p\ 2 ~^ R- e [^a ] 2 + ie, in agreement with 
Eq. (I2.2p . Instead, in giG(s) ao p, with P one of the lightest pseudoscalars, 
M ao should appear with a negative imaginary part to guarantee that, in the 
zero-width limit, the sign of the imaginary part is the same as dictated by 
the — ie prescription for masses squared of the intermediate states. Such 
analytical extrapolations in the masses of external particles are discussed in 
Refs. [501452] . The same applies of course to the case of the /o(980) resonance. 
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3 Results 



In this section we show the results that follow by applying eq. (12. lip to the 
different channels characterized by the quantum numbers IG, as given in the 
list from eq. ( 12. 6ft to eq. (12. 10ft . As discussed after eq. ( 12. 12ft we consider 
values for the subtraction constant such that they are negative and not very 
large in modulus (|ai| < 3). In this way, the resonances generated might be 
qualified as dynamically generated due to the iteration of the unitarity loops. 
The pole positions and couplings of the /o(980) and ao(980) resonances are 
given in Tables [3] and HI respectively, and they correspond to those obtained 
in the meson-meson S-wave amplitudes used, see appendix [A] We present 
the results for each of the channels with definite IG separately. 




1200 1400 1600 1800 

y/s (MeV) 



0.6 




0.0 1 1 1 1 1 1 1 

1200 1400 1600 1800 

y/s (MeV) 



Figure 3: Modulus squared of the foK — > f K (left) and a§K — > aoK (right) 
S-wave amplitudes for ai = —0.5 . The points correspond to the energies 
where the amplitudes have been actually calculated. 



10 



3.1 7 = 1/2 



First we show the results for the 7 = 1/2 sector that couples together the 
channels /o(980)if and a (980)K. We show the modulus squared of the 
foK — > f K and a^K — > a K S-wave amplitudes in the left and right panel 
of Fig. [31 respectively. We obtain a clear resonant peak with its maximum at 
1460 MeV for a-^ around —0.5, that corresponds to the nominal mass of the 
i^(1460) resonance [14j. The results are not very sensitive to the actual value 
of ai but the position of the peak displaces to lower values for decreasing ai 
and the width somewhat increases. The visual width of the peak is around 
100 MeV, although it appears wider in a^K — > a^K scattering. In refs. [53154"] 
a larger width of around 250 MeV is referred. One has to take into account 
that the channel K*(892)n is not included and it seems to couple strongly 
with the i^(1460) resonance [H]. It is also clear from the figure that the peak 
is asymmetric due to the opening of the foK and a^K thresholds involved. 
Taking into account the relative sizes of the peaks in the left and right panels 
of Fig. [3] one infers that the -fT(1460) couples more strongly to JqK than to 
a K, with the ratio of couplings as \gf K / 9a K\ — (^§y) 1//4 — 1-4 . 

3.2 1=1 

We now consider the 1 = 1 case. As commented in the introduction two 
broad resonances are referred in the PDG, the 7r(1300) and 7i~(1800). In our 
amplitudes we find quite independently of the value of ai that the 00(980)77' 
channel is almost elastic. This is due to the fact that the interaction kernels 
T(ao?7 / — > a rj) and T(aofj' — > foil) are much smaller than the rest of kernels, 
typically by an order of magnitude. This happens because the kernels are 
dominated by the threshold region. However, the threshold for ao (980)77' is 
much higher than the thresholds for the other two channels. In this way, 
for the inelastic processes involving the a i]' channel, even at threshold for 
one of the channels, there is always a large three-momentum for the other 
channel and the kernel is suppressed. Of course, this does not apply for the 
0077' elastic case where the kernel has a standard size and produces around 
1.8 GeV a strong resonant signal that could be associated with the 7r(1800) 
resonance. To reproduce the mass value given in the PDG [H] for this 
resonance, 1816 ± 14 MeV one takes a\ for ao7/ around -1.3. The visual width 
of the peak is around 200 MeV, close to the width quoted in the PDG [H] of 
208 ±12 MeV. The other two channels couple quite strongly between each 
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1000 1200 1400 1600 1800 2000 1000 1200 1400 1600 1800 



y/s (MeV) v/i (MeV) 

Figure 4: Modulus squared of the dor]' — >■ a r/ (left) and ao?7 —> aoi] (right) 
S-wave amplitudes. For the former a\ = —1.3 and for the latter a% = —2.0, 
see the text for details. The notation is as in Fig. [3j 

other and typically give rise to an enhancement between 1.2-1.4 GeV when 
varying a% equal for each of them, which could be associated with the 7r(1300). 
However, for |ai| between 1 and 1.8 a too strong signal in the aof] threshold 
originates. For |ai| below 1 the resonant peak in the \T(aor] — > a i])\ 2 lies 
around 1.4-1.5 GeV, somewhat too high for the 7r(1300) resonance [H]. This 
is why we show in Fig. H] the modulus squared of a ?7 —> a Q r] for a x = —2 
where a peak close to 1.2 GeV is seen with a width of around 200 MeV. One 
can also see the strong effect of the a rj threshold at around 1.52 GeV. Its 
size is rather sensitive to the actual vale of |ai| when this lies between 1 and 
1.8. There is the interesting fact, which is independent of the value of a 1; 
that there is no signal for 7r(1800) in the aorj system nor signal of the peak 
at 1.2 GeV in the dorf. We have also checked that this is also the case for 
the foil state, that is, it does not couple with the 7r(1800). This is another 
reflection of the fact that the dorf tends to decouple from the other states. 
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Figure 5: Modulus squared of the /or/ — > forf (left) and /or/ — > foV 
S-wave amplitudes. For the former ai = —1.25 and for the latter a\ = 
see the text for details. The notation is as in Fig. [3j 



1800 



[right) 
-0.8, 



3.3 1 = 

We move next to the 1 = system where the for], a^ir and /or/ couple. Here 
occurs similarly to I = 1, so that the much higher for]' channel mostly decou- 
ples from the other two channels. We then proceed similarly and distinguish 
between the subtraction constant a\ attached to a i]' and to the other two 
channels a^it and /o 7 ?- F° r a i around —1.2 one obtains a resonance of the 
dorf channel at a mass of 1835 MeV, in agreement with that quoted in the 
PDG for the X(1835), 1833. 7±6.1 ±2.7 MeV. This is shown in the left panel 
of Fig. [5] where the modulus squared of the /o(980)r/ —> /o(980)r/' S-wave 
amplitude is shown. The width of the peak at half its maximum value is 
around 70 MeV, in good agreement with the width given in the PDG for the 
X(1835) of 67.7 ± 20.3 ± 7.7 MeV. We consider next the other two coupled 
channels, ao(980)7r and /o(980)r/. We obtain a clear resonant signal with 
mass around 1.45 GeV for |ai| < 1. This is shown in the right panel of 
Fig. where the modulus squared of the / (980)r/ — > / (980)?7 is given for 
a\ = —0.8. It is not possible to increase further the mass of this peak by 
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varying ai. An important fact of this resonance is that it does not couple 
to the a 7i channel. E.g. the analogous curve for the modulus squared of 
the ao(980)7r — >■ a (980)7r S-wave in the 1.4 GeV region is absolutely flat. By 
considering the inelastic process /o(980)7/ — > a (980)7r we estimate a coupling 
to the latter channel more than 14 times smaller than to /o(980)r/. Because 
the 7/(1405) resonance couples mostly to a (980)7r [H] we then conclude that 
the generated resonant signal around 1.45 GeV should correspond to the 
77(1475) . Its form is rather asymmetric due to the opening of the /o(980)7/ 
threshold, with a width at half the maximum of its peak of around 150 MeV. 
The width quoted in the PDG [14J is 85 ± 9 MeV. It is also known that the 
77(1475) couples strongly to K*(892)K + c.c, a channel not included in our 
study. The threshold for this channel, at around 1.39 GeV at the decreas- 
ing slop of our present signal, should certainly modify its shape. For higher 
values of |ai| the peak tends to become too light in mass compared with the 
77(1475). For the ao(980)7r — > a (980)7r reaction one also appreciates a strong 
a (980)7r threshold effect at around 1.16 GeV. No resonance around the mass 
of the ^(1295) is observed. 

0.40 1 




0.00 1 1 

1000 1200 1400 1600 1800 2000 
^/s (MeV) 



Figure 6: Modulus squared of the 1 = 3/2 a K — > a K S-wave amplitude 
with a\ = —0.5. The notation is as in Fig. |3j 
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Resonance 




Width (MeV) 


Comments 


#(1460) 


i = h 


r>ioo 


\gf K/g ao K\ ^ 1.4 


tt(1800) 


i g = i- 


r ~ 200 


a r]' elastic 


tt(1300) 


i g = i- 


r>2oo 


^o 71 "; foV coupled channels 


X(1835) 


I G = + 


r ~ 70 


forf elastic 


7?(1475) 


I G = 0+ 


r ~ 150 


forj elastic 


Exotic 


1 2 


r ~ 200 


a K threshold 



Table 1: Resonances resulting from our study. For more details see the discussions 
of the results in the text. 

3.4 Exotic channels 

Regarding the exotic channel with I = 3/2 we find an interesting result. 
Our amplitude gives rise to a clear resonant structure at around 1.4 GeV for 
| cfci | < 1.5. We show the modulus squared of T{a^K ± — > a^K^), because 
the clqK ± states are purely I = 3/2, for a\ = —0.5 (the same value used 
before in Fig. [3] studying the 1=1/2 case) in Fig. [6j One also observes that 
the shape of the resonance peak is asymmetric with a clear impact of the 
clqK threshold. Our results for |ai| < 1 tends to confirm the predictions of 
Longacre [55] that studied the KKn and KKK system and concluded that 
the exotic I = 3/2 J p = 0" KKK system was resonant around its threshold 
due to the successive interactions between a K, K and a ir. For \ai\ > 1 we 
find that the resonance shape in \T{a^K ± — > a^K^l 2 progressively distorts 
becoming lighter and flatter. Let us notice also that the a$K system was 
not isolated in the two experiments quoted in the PDG where the 7 = 1/2 
#(1460) was observed ^SB. 

The other exotic channel with 1 = 1 and G = +1 involves the isovector 
a iT state. Whether a resonance behavior stems at around 1.4 GeV depends 
on the actual value of a%. For |ai| < 1 the enhancement near 1.4 GeV is much 
weaker and is overcome by the cusp effect at the aovr threshold. For larger 
values of |ai| the resonant signal is much more prominent. No such resonance 
has been found experimentally, e.g. in peripheral hadron production [56] . so 
that | ai| < 1 should be finally taken. 

In Tabled] we collect all the resonances found in our study for the different 
quantum numbers discussed. 
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4 Summary and conclusions 



In summary, we have presented a study of the S-wave interactions between 
the scalar resonances /o(980) and ao(980) with the lightest pseudoscalars (tt, 
K, 7] and rj') in the region between 1 and 2 GeV. The different channels stud- 
ied comprise those alike the rj, K and tt, and the exotic ones with isospin 3/2 
and 1, the latter having positive G-parity. First, interaction kernels have been 
derived by considering the interactions of the external pseudoscalars involved 
in the reaction with those making the scalar resonance. We take advantage 
here of previous studies that establish the / (980) and a (980) as dynami- 
cally generated from the interactions of two pseudoscalars, so that no free 
parameters are introduced in their calculation. Afterwards, the final S-wave 
amplitudes are determined by employing techniques borrowed from Unitary 
Chiral Perturbation Theory. Interestingly, we have obtained resonant peaks 
that for the non-exotic channels could be associated with the pseudoscalar 
resonances K(U60), vr(1300), vr(1800), r?(1475) and X(1835), following the 
notation of the particle data group. The resonances that come out from this 
study can be qualified as dynamically generated from the interactions be- 
tween the scalar resonances and the pseudoscalar mesons. This establishes 
that an important contribution to the physical signal of the resonances just 
mentioned has a dynamical origin. The exotic 7 = 3/2 channel could also ex- 
hibit a resonant structure around the a^K threshold, in agreement with the 
behavior predicted by Longacre [55] twenty years ago. However, larger values 
for the subtraction constant |ai| tends to destroy this resonant behavior. No 
signal of the intriguing i] (1405) resonance is obtained. 

This approach should be pursued further by including simultaneously to 
the interaction between the scalar resonances and the pseudoscalar mesons, 
considered here, those arising from the lightest vector resonances with the 
same pseudoscalars in P-wave. In this way, both pseudoscalar and axial 
resonances will be studied together. 
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A Meson— meson unitarized amplitudes 

We use the N/D method fTS] to unitarize the different isospin channels am- 
plitudes for meson-meson scattering, which are fitted to data, and then used 
in the vertex of the triangle loop. From these amplitudes, once fitted, the po- 
sition of the poles can be found (we use here the /o(980) and ao(980), but we 
also check for the appearance of the other scalars, a and k). As mentioned, 
the amplitudes are unitarized through 

T I = (1 + V I -G)- 1 -V I , (A.l) 

which is analogous to eq. (12.111) but now for the pseudoscalar-pseudoscalar 
scattering. The symmetric matrix Vi (the analogous one to Tig m eq. (12.111) ) 
collects the S-wave pseudoscalar-pseudoscalar tree-level amplitudes obtained 
from the lowest order Chiral Lagrangians including resonances as well. The 
matrix G is a diagonal matrix that contains the meson-meson loop prop- 
agator (the same expression as given in eq. (I2.12p can be used with the 
appropriate replacement for the masses involved.) 

The lowest order chiral Lagrangian at leading order in large N c which 
also includes the 771 is [571459] 

L 2 = f l{d^3^U) + f -^(x ] U + X U ] ) - l -Mlr,l , (A.2) 

U{4>) = exp(iV2$/U), $ = E7/f^ Ao = V3 l3 ' ( A - 3 ) 

and Aj, i — 1, . . . , 8 the Gell-Mann matrices. The 77s and rji fields mix to give 
the physical rj and 7/ 

( 1 _ ( cos ^ sin^ \ / 771 \ 
I 77 J \ — sin 9 cos 9 J \ 7/ 8 J 

The mixing angle 9 is taken as sin 9 ~ — |, 9 ~ —20° [60J. 
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Figure 7: Eq. (I A . 5 P in terms of Feynman Diagrams 



In the same spirit as in ref. [13] the explicit exchange of J PC = ++ 
scalar resonances is incorporated and calculated from the leading order chiral 
Lagrangians of ref. |61j. The appropriate Lagrangians are 

(A.4) 



Ls 8 


= c d (Sgu^) + 


Cm (S 8 X+) 








= c d S% (u^) + 


CmSl (X+) 






X+ 


= u'x u + ux^u 








U{x) 


= u(x) 2 , Up = 


iu^dpUv) = 










/ 0,0 , fs 




K* + 




s 8 






K*° 


K*° 








V *r 




> 



and Si is a scalar SU(3) singlet. The interaction kernels obtained from La- 
grangians (1A.2|) . (IA.4[) and (1A.5|) can thus be written as 



V: 



Vii 
(R) 



V {c) + V {R) 
cacti 



(A.5) 



'.i 



M 2 R -s 



where C means contact term and R resonance exchange. This is represented 
diagrammatically in Fig. [7J In what follows, we give explicit formulae for 
the contact kernels from the chiral Lagrangians in eq. ( 1A.3I) for the different 
isospin channels. We also give the couplings a{ for each scalar resonance. For 
/ = we include the superscript (8) or (1) in «j to distinguish between the 
octet and singlet contributions, respectively. For the rest of isospins there is 
no singlet contribution. 
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/ = 



V (C) 



2s - ml 

s 



V {c) - 

= -— ^ 

V KK^KK 4 



(c) 2 3s - 2m^ - 

(c) = 2m 2 K + ml (A6) 



a 



(8) 



a 



(8) = 
KK 

a {8) = 

tfri 



c d s + 2(c m - c d )ml 
f 2 

c d s + 2(c m - c rf )m|- 
8c m (m 2 K - m 2 ) 



(A.7) 



(i) 



a 



(i) 



= 

Tfri 



r- cgs + 2(c m - c d )ml 

V f 2 

c d s + 2(c m - c d )m|- 
2v/2/ 2 

^ 3c d (s - 2m 2 ) + 2c m {2m 2 K + ml) 



(A.8) 
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a 



KT} 



®KK 



2 c d s - c d (ml + m 2 ) + 2c m m\ 

"7! ft 
c d s - 2(c d - c m )m\ 

f 2 

2 c d s - c d {m\ + ml,) + 2c m m\ 



f 2 

J 7T 



(A.12) 



/ = 3/2 



^Ktt 



in 



K 



mz 



V; 



(A.13) 
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Figure 8: Experimental data and our fits, as explained in the text. 

With the amplitudes calculated for meson-meson scattering, we per- 
form several fits, e.g. by changing the value of the highest y/s fitted from 
1.2 to 1.4 GeV and by imposing that several subtraction constants for the 
pseudoscalar-pseudoscalar channels are equal, so that we can calculate the 
pseudoscalar-scalar kernels with different inputs, and then check the inde- 
pendence of our results. We only show our main fit since all the other fits 
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that we obtained give rise to similar results that would not change our con- 
clusions. In this fit the highest value of y/s considered is 1.4 GeV. For the 
octet of scalar resonances we take the values of the parameters q, c m and 
M 8 from ref. (60], where M 8 , the mass of this octet, is around 1.3-1.4 GeV. 
For definiteness, q = c m = 22.8 MeV and M 8 = 1.4 GeV. The parameters 
for the singlet resonance exchange, q, c m and Mi are left free, with the 
latter the mass of the singlet scalar resonance. Regarding the subtraction 
constants in the unitarity loop function of the different channels [13] (they 
play the analogous role of a\ in eq. (I2.12p but for pseudoscalar-pseudoscalar 
scattering), we take the most general situation compatible with isospin sym- 
metry. Adopting the same argument as in the appendix A of ref. [62] from 
SU(3) to SU(2), the subtraction constants corresponding to the same pair of 
pseudoscalars should be the same in the different isospin. In this way, the 
subtraction constant for KK both in I — and / = 1 is taken with the 
same value. On the other hand, for a given isospin, we also put constraints 
on the subtraction constants associated with non-relevant channels. In this 
way, the Kt] subtraction constant in I = 1/2 is kept equal to that of Ktt 
and, similarly, for 1=1 the tti]' subtraction constant is put equal to that 



of Of course, we have checked that smoothing these constraints does 

not affect the results of the fit. In this way, we finally have six independent 
subtraction constants for tttt, KK, rjr], Kir, Krj' and ttt]. There is also a nor- 
malization constant for the data on an unnormalized ttt] event distribution 
around the ao(980) resonance that is required for each fit. The results of the 
fit compared to experimental data are shown by the solid line in Fig. [8] and 
the values of the fitted parameters are given in Table [2j 

The set of experimental data included in the fits for 1 = comprises 
the elastic tttt phase shifts, 5^^^, from refs. [SBTfTD] . the phase shift for 
tttt KK, 5nK-*KKi an d (1 — V 2 )/^ from refs. [7Tj|72], where rj is the elastic 
parameter for the tttt — > tttt I = S-wave. With respect to / = 1/2 we fit the 
elastic ttK phase shifts, 5ktt^Ktt, from refs. [73TTT6] . Finally, we include an 
event distribution of ttt) around the ao(980) resonance mass from the central 
production of ttttt), ref. [77], fitted like in ref. [TT1IT2] . 

Once the fits are performed we look for the poles of the scalar resonances 
cr, k, / (980) and a (980) in the unphysical Riemann sheets continuously 
connected with the physical one. Notice that only the /o(980) and a (980) 



^ tl We have checked that the 7rr/ channel tends to decouple of the ttt] and KK channels 
in I = 1 in the region of the ao(980). 
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Parameter 


Value 


c d (MeV) 


18 ± 1 


cm (MeV) 


23 ±4 


Mi (MeV) 


1100 ±20 




-0.98 ±0.10 


a KK 


-1.00 ±0.20 


CLrjrj 


+0.04 ±0.22 




+0.17±0.05 


a,Kr)' 


-3.53 ±0.13 




-2.55 ±0.37 



Table 2: Fitted parameters for the main fit. The x 2 /d.o.f. is 0.96. The fits 
are obtained employing the program MINUIT [S3] . 



Resonance Re^/s (MeV) hny/s (MeV) 



a 


466 


235 


K 


698 


294 


/o(980) 


987 


18 


a (980) 


1019 


33 



Table 3: The pole positions of the scalar resonances obtained from the main 
fit are given. 



Resonance g K g (GeV) \Qkk (GeV)| 



Jo (980) 


-3.72 4 


- 1.18i 


3.90 


a (980) 


-4.11 4 


- 1.59i 


4.41 



Table 4: Couplings of the /o(980) and ao(980) resonances to KK (with 
definite isospin). These couplings are calculated from the residues of the 
corresponding pole, see e.g. [T3~|l6~4"] . 
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poles are actually required for evaluating the pseudoscalar-scalar scattering 
kernels in section |5J The a and k are given for completeness. They are 
related to the /o(980) and ao(980) resonances giving rise to a nonet of light 
scalar resonances jH]. Other poles around 1.4 GeV also appear that we do 
not include here. The pole positions are given in Table [3j The couplings of 
the /q(980) and ao(980) to KK, used in this work, are collected in Table SJ 



B S-wave projection of C3 and C±(Ml) 

The three- and four-point Green functions C3 and C^mf) are defined in 
eq. (12. 5p . Here, we consider the more general case with arbitrary internal 
masses and from the very beginning the S-wave projection is worked out. 
Both functions are finite. 




Figure 9: Feynman diagrams for C3 (left) and C^{m^) (right). 



We first consider C3, left diagram of Fig. [9J and follow its notation with 
t — (pi + P2) 2 (note that all four-momenta are in-going). We also introduce 
two Feynman parameters U\ and u 2 and the relative angle 9 between the 
initial and final pseudoscalars, so that 

Cs= 2j_ 1 dC ° s9 J (2tt) 4 ((£ + Pl y - m\){P - ml)((£ + Pl +p 2 y - m§) 

= oo - 2 / / du 2 I d cos 9 \p\u\ + p\u 2 2 + 2pip 2 u\u 2 
<j27T J q J q J_ 1 

+ {m\ - ml - p\)ui + (ml - ml - p\ - 2p 1 p 2 )u 2 + ml + ie] 1 . (B.l) 
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Then, cos# is introduced by taking into account that p\p 2 — (t — p\ — pi)/ 2 
with t = g 2 — |p| 2 — |p'| 2 + 2|p||p'| cos#, with q° = p\ — the difference 
of energies between the initial and final scalar resonances. We perform the 
angular integration and introduce the parameter £ 2 as u 2 = u\^ 2) so that 

° 3 = 1 1 I ' f |log(1 + *> ~ log( - 1 + *>> ' (R2) 

where 

; ^ T 2 2 1 2 2 1 / 2 I 1 2 I ' 1 2 2 2\ 

^ = i I, nn x— [P2 u i +PiU 2 + Uiu 2 {q - p - P -P1-P2) 

2|p||p 1(1 - U\)u 2 

1 ^2 2 2\ 1 / 2 2,2 2,||2,| / 12\ , 2 ■ 

+Mi(m 1 - m 3 — p 2 J + -U 2 (^2 _ m l + £>2 _ % + |Pl + IP I ) + m 3 _ l£ \ 

(B.3) 

For the four-point function C^m 2 ), right diagram of Fig. |9l one has 
^ , 2, * f +1 r a f d 4 £ 1 

G4 m- = - GSCOStf / — ~7T~ tt; ON . — =7 

' (B.4) 



((£ -Pa -Pi) 2 ~ ml)((£-p 4 ) 2 - ml) 

In this case there is no ambiguity if instead of performing the cos 9 integration 
one directly calculates the related integration over t = {P3+P4) 2 by taking into 
account that dt = 2|p||p'|dcos0 (ambiguities could arise for s such that the 
product IpIIp'I becomes complex. The particular integration to be performed 
here is not affected by such problem, see below.) We also introduce three 
Feynman parameters ui, u 2 and u 3 so that 

1 rl ru\ i-u 2 rt + 

C 4 (m 2 ) = / dui / du 2 I du 3 / dt [p\{l -u 2 + u 3 ){u 2 - u 3 ) 

647r 2 |p||p / | J J J J t 

+2p 1 p 3 (u 1 - u 2 ){u 2 - u 3 ) + 2pip 4 (l - u 2 )(u 2 - u 3 ) - m\u 3 

+pl{l -Ui + u 2 )(ui - u 2 ) + 2p 3 p 4 (u! - u 2 )u 2 +pl(l - u 2 )u 2 

-m\{u 2 - us) - m 2 (wi - u 2 ) - m\{l - u\) - is] 1 . (B.5) 

In terms of the variable t the previous integral is of the from J dtjiat + 6) 2 
so that the t-integration can be done straightforwardly without problems 
in its analytical extrapolation. The resulting M 3 -integration is of the form 
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J du^f [us(ul+ Pus+j)} that can also be done straightforwardly by factorizing 
the second-order polynomial in the denominator. Our final expression for 
£4(7774) is 




L 



o 1-6 



[V>(f + ,u 2 )-V(t+,o) 



^(t_,«2)+V(t-,0)] 



(B.6) 



where 



«3 + P U 3 + 7 
P?(u3 + /9«3+7) 



V>(*,us) 



lQg^3 log(M 3 ~ 2/l) _ log(M 3 ~ j/g) 

2/12/2 2/1(2/1-2/2) 2/2(2/1-2/2) 

(«3 - 2/l) («3 -2/2) , 

7774(1 - Mi) + m^-ui - -u 2 ) +Pa(-1 + - U2) + mju2 



+ u 2 (s(—l + ui) + ^2(^2 - Mi)) - (ml -ml+ pl)u 3 + plul 
+ (s(l - Mi) + (pi +pl~ t)u! - (pi +pj- t)u 2 )u 3 - is , 



The S-wave projection of the three- and four-point functions C 3 and 
(74(777.4) was also obtained for some kinematical regions and values of 7774 
from ref. [78], wherever the latter could be applied. In such cases our results 
and ref. [78] agree. 
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